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Abstract 

I present the junction conditions for F(R) theories of gravity and their 
implications: the generalized Israel conditions and equations. These junc- 
tion conditions are necessary to construct global models of stars, galaxies, 
etc., where a vacuum region surrounds a finite body in equilibrium, as 
well as to describe shells of matter and braneworlds, and they are stricter 
than in General Relativity in both cases. For the latter case, I obtain the 
field equations for the energy- momentum tensor on the shell/brane, and 
they turn out to be, remarkably, the same as in General Relativity. An 
exceptional case for quadratic F(R), previously overlooked in the litera- 
ture, is shown to arise allowing for a discontinuous R, and leading to an 
energy-momentum content on the shell with unexpected properties, such 
as non-vanishing components normal to the shell and a new term resem- 
bling classical dipole distributions. For the former case, they do not only 
require the agreement of the first and second fundamental forms on both 
sides of the matching hypersurface, but also that the scalar curvature R 
and its first derivative VR agree there too. I argue that, as a consequence, 
matched solutions in General Relativity are not solutions of F(i?)-models 
generically. Several relevant examples are analyzed. 

PACS: 04.50.Kd; 04.40.Dg 
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1 Introduction 



Alternative theories of gravity, such as models based on an F(R) Lagrangian, 
have received considerable attention in the past years, see e.g. [TJ [2j [3]. The 
basic fundamental solutions of standard General Relativity (GR), such as the 
Schwarzschild and Kottler spherically symmetric exteriors, or the Friedmann- 
Lemaitre-Robertson- Walker cosmological models, are also solutions of the F(R)- 
theories. In the latter case, however, the energy-momentum content of the space- 
time differs from that of GR. Concerning the former case, and as an illustrative 
example of the greater richness of F(R) theories, the Birkhoff- Jensen theorem 
on the uniqueness of the exterior (vacuum) solution of a spherically symmetric 
gravitating system no longer holds and thus there arise many other (even non- 
static) vacuum solutions different from the Schwarzschild solution [2]. Similarly, 
many possible interior solutions (be they perfect fluids or more general matter) 
are feasible in the extended theories. It is tacitly understood, in general, that 
the GR solutions are a very good approximation to the F(R) solutions at least 
in some regimes such as low energy ones [2j. 

In this paper, the question of whether or not global models in GR describing 
(say) stars in equilibrium or collapsing bodies that form black holes — in both 
cases all the way from the centre up to infinity — remain as viable valid models in 
the generalized F(R) theories will be addressed. The related problem of inserting 
cavities into standard cosmological models, which requires a similar analysis, has 
been recently considered in jl] . The conclusions obtained herein agree with theirs. 
In both cases, one needs the proper junction conditions through the hypersur- 
face separating the exterior and the interior of the global gravitational field, in 
particular, through the world-surface of finite bodies such as stars. These junc- 
tion conditions have been obtained in [5] via a calculation performed in Gaussian 
coordinates relative to the matching hypersurface, and were later used in jl]. 
The calculation in [5] allows for cases with branes or thin shells at the matching 
hypersurface. However, the field equations for the energy-momentum tensor on 
the brane (formulas (jl]) and (JSJ) below) were not derived in j5], and they do not 
seem to have been written down hitherto. Remarkably, they are actually identical 
with those of GR, even though the derivation as well as the expression for the 
energy- momentum tensor are quite different in F(R) theories. This leads to some 
speculations about its possible universal character. 

To avoid loading the paper with heavy calculations, I derive the correct junc- 
tion conditions for F(R) gravity in an Appendix by using tensor distributions 
[5] [7]. The main results are succinctly presented in section [2] For non-linear 
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functions F(R), the junction conditions always require continuity of the trace 
of the second fundamental form of the matching hypersurface as well as, geneti- 
cally, the continuity of the scalar curvature R. It turns out that matter shells and 
braneworlds cannot be umbilical hyper surf aces, and the brane tension is propor- 
tional to the discontinuity of the derivative of R, in contrast with the GR result, 
see however (HI E] and references therein. 

Very surprisingly, there also arises an exceptional case which seems to have 
been overlooked in the literature and where a discontinuous scalar curvature R 
is allowed. This exceptional case can only arise in theories with a quadratic 
function F(R), and leads to quite a different matter content on the brane or 
thin shell so that, for instance, components normal to the matching hypersurface 
can arise as well as a completely unexpected new term, described by (jUJ) and 
fl39|) below, which resembles those of classical "dipole distributions" — a kind of 
Dirac-delta-prime distribution. The interpretation of such a new term is quite 
unclear. This is considered in subsection 12.1.21 where the general field equations 
for the energy-momentum content of the shell in this exceptional situation are 
presented. 

The case without braneworlds or thin shells, so that the curvature tensor 
distributions do not possess singular parts, is then considered in subsection 12.21 
and their implications in terms of the energy-momentum tensor quantities - 
equation (fl3|) below — are derived. This seems to be also new and again adopts 
the very same form as in GR, despite this not being trivial at all. This allows 
to prove, for example, the important result that the matching hypersurface for a 
compact perfect fluid is always defined by the vanishing of the pressure, and then 
the fluid has also vanishing energy density there or at least it becomes tangent 
to the matching hypersurface — see subsection 12.2.11 

The junction conditions turn out to be more restrictive that in the GR case, 
as they will impose — in addition to the same conditions as in GR — differen- 
tiability of the scalar curvature R across the matching hypersurface. This will 
lead to a simple but general proof that, in general, GR matched solutions will 
not be solutions of generalized F(R) gravity theories, as analyzed in section [31 
In particular, for example, the Oppenheimer-Snyder collapsing star [S] to form a 
black hole is not a solution of any F(i?)-theory. The corresponding complemen- 
tary matching in the sense of [10J, which describes the Einstein- Straus vacuole 
[IT], see [12] . is also impossible for non-linear F(R) as recently demonstrated in 
[I]. The constant-density interior Schwarzschild solution matched to its exterior 
is not a solution either, nor are the vast majority of the multiple static spherically 
symmetric solutions with a perfect fluid interior matched to the exterior vacuum 
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Schwarzschild solution — see also, in this respect, the discussion in [T31 E]. I 
also consider more general cases with dynamical and radiating exteriors, and the 
conclusion is the same: GR matched solutions are not solutions of the extended 
F(R) theories generically. A very particular model that is a global solution of 
both GR and F(R) theories is found, however, given by a Robert son- Walker in- 
terior with radiation equation of state matched to the Vaidya radiating exterior 
solution [13]. The complementary matching gives rise to a radiating vacuole in 
an expanding universe. 

I must mention that there have been several papers discussing the different 
problems that arise concerning exact solutions of F(R) theories, such as [TBI [TT] 
[T8"] and many references therein, where curvature singularities appear or there 
arise impediments to have strong gravitational fields in spherically symmetric 
stars. The difficulties shown herein are, nevertheless, of an intrinsic different 
nature as I am just concerned with the question of whether or not matched GR 
solutions can be solutions of the extended theories. 

2 The matching conditions 

The appropriate framework to study the matching of two different spacetimes 
across a timelike hypersurfac^j] is that of tensor distributions [BJ [7] , because the 
proper junction conditions follow from analyzing the singular parts of some cur- 
vature tensor distributions and, for the proper case with no brane/shell, by de- 
manding that they vanish. This leads in GR to the standard Darmois and Israel 
matching conditions [20l |2T| requiring the agreement, on the matching hyper- 
surface, of the first fundamental form, while supplying a formula (respectively, 
field equations) for the brane/shell energy-momentum content in terms of the 
discontinuities of the second fundamental forms (resp. of the energy-momentum 
tensors) inherited from both sides of the space-time . In particular, there is no 
shell/brane when the second fundamental forms agree, in which case the junction 
conditions are equivalent, in a certain sense, to the Lichnerowicz C 1 conditions 
on the metric components in admissible coordinates [22], [20], |2TJ . For a summary 
of the junction conditions and standard references see section 3.8 in [23] and [19] . 
I present a derivation of the junction conditions following this general method in 

lr The whole study can be performed for null hypersurfaces, and in general for hypersurfaces 
changing their causal character from point to point |19j . However, to keep the presentation 
as simple as possible and make the main point plain I have preferred to restrict myself to the 
important case of a timelike matching hypersurface. 
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the Appendix, and I only give in this section a summary of the results, stressing 
the new formulas that were not found before and the appearance of a particular, 
extraordinary case, for theories with a quadratic function F(R). 

Following the notation used in the Appendix, let £ be the timelike matching 
hypersurface, the unit normal to S, h^ u = g^ u — n^n u its first fundamental 
form and K^ u its second fundamental form. This last object, as well as others, 
may have a jump across S. For any function / its jump across £ is denoted by 

[/]■ 

2.1 Case allowing for matter shells or branes 

It is shown in the Appendix that in general F(R) theories the following condition 
must hold 

[K%] = . (1) 

Therefore, the trace of the second fundamental form must always be continuous 
across S, even for cases with shells of matter or branes. This is quite different 
from the GR case, and forbids the use of umbilical hypersurfaces — characterized 
by = fh^y — one of the most common cases to describe braneworlds in GR. 

Now the analysis splits into two possibilities, depending on whether or not 
F"'(R) =0. 

2.1.1 The generic case: F"'(R) ^ 

In this CcLSG, clS shown rigorously in the Appendix, the requirement 

[R] = (2) 

is unavoidable. Hence, the scalar curvature must always be continuous across E, 
even for cases with shells of matter or branes. 

The energy-momentum tensor in the brane or thins shell reads (see Appendix) 

«V = -F'(R S ) [K^] + F"(R s )nP [V P R] V, ^ V = 0. (3) 

Thus, the following result holds: 

The proper junction conditions allowing for shells of matter or 
branes in F(R) theories with F"'(R) ^ are the agreement of the 
first fundamental forms on both sides of the matching hypersurface 
together with (T5]) and ([1]). The energy- momentum content of the shell 
or brane is then given by formula (j3J). 
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It should be noted that the first summand in fl3]) is traceless due to ([T]), so that 
the trace of the energy-momentum singular part (sometimes called the "brane 
tension") reads simply 

KT% = (n-l)F"(Rx)n p [V P R] 

and is fully determined by the discontinuity of the normal derivative of R across 
E. This is in sharp contrast with the GR case, where the brane tension is the 
discontinuity of K p p , which is always continuous now. In particular, for S with 
only brane tension A so that kt^ = \h^ v one must have [K^] = and the tension 
is A = F"(Rz) n p [V P R]. 

Concerning the equations satisfied by the energy-momentum tensor in the 
brane, one has from the Appendix 

(K+ + tf - )t" = 2nV [T Plt ] (4) 
V%„ = -n p h\ [T pa ] . (5) 

As far as I am aware, these relations are new and, most remarkably, identical 
with the GR case — corresponding to ( 1271) and (125]) via the Einstein equations. 
It arises the idea that they may be universally valid for diffeomorphism-invariant 
theories , but I do not know of any result along these lines. 

2.1.2 The exceptional case: F"'(R) = 

As proven in the Appendix, this case allows for a discontinuous R, so that [R] ^ 
is possible and the energy-momentum tensor distribution acquires a singular part 
with two terms 

= T+9 + T~(l - 0) + (r pv + T p n u + r u n p + rn p n v ) 5 s + 

where 

«V = -{1 + a(R + + R-)}[K pu ] + a {2ah pu - [R](K+ V + K~ u )} , n%, = 0, (6) 

KT ll = -2a%[R], n% = 0, (7) 
kt = 2a[R}K%, (8) 

k!^, = 2^- (9) 

Here 

a = l ~F"{R) 
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is a constant, a is a function on E defined in (1551) and the delta-prime-type 
distribution Q is defined in ( 1391) . 
Thus, we now have: 

The proper junction conditions allowing for shells of matter or 
branes in F(R) theories with F"'(R) = are the agreement of the 
first fundamental forms on both sides of the matching hypersurface 
together with (JT|). A discontinuous R is permitted and the energy- 
momentum content of the shell or brane is given by formulas (J6H9J) . 

The discontinuity of the energy-momentum tensor can be computed from ( l3Tj) 
with f )34|) . leading easily to 

= (1 + 2aR^)[G„ v \ + 2a{[R]R% + g pu [V p V p R] - [V P V U R]} 

which on using f )36|) . (12?]) . f )28|) and fl25|) provides in this case 

Kn p n v [T pv \ = -(1 + 2aR- s )K^[K pa ] + 2a { [R]R* u n p n u + aK p p + W P W P [R] } , 
Kn v h\\T pv \ = (l + 2aR s )V P [K pfl }_ _ 

+ 2a { \R\n v h p p R pv - V p a + [K pp \ V P Rx + K^V P [R] } . 

Combining these with (EHH]) and ( )29ll30l) one obtains 

n u h%[T p „] = -V^v-ifVe-V. 
nV[^]- v C + V\ = 2-[R}(R*nW + K^K s ^) 

K 

= -[R]{R^-TZ+(K%) 2 + K^K-^). 

It may be observed that the case a = is simply GR. If on the other hand one 
sets [R] = then the generic case is recovered under the assumption fl34l) . 

The appearance of the last term (Q is most remarkable, and very surprising. 
The classical interpretation of such terms (for instance in electromagnetism) is 
that they describe a distribution of dipoles on E. This is usually thought to 
be physically not viable in Gravitation due to its attractive character and the 
positivity of masses. The interpretation of such new terms in this exceptional 
case is thus open and may give rise to new possibilities for describing quite exotic 
braneworld scenarios. 
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2.2 Case without shells or branes 



Consider now the proper matching case where no shells of matter or branes are 
allowed, so that T^ u can have, at most, jump discontinuities. From the Appendix 
we know that in this situation (T2J) must hold and also 

[K„ v ] = 0, n" [V P R] = 

thus 

[V M i?] = 0. (10) 

This immediately informs us that the GR junction conditions are necessary, that 
is to say, the first and second fundamental forms must agree on E from both sides, 
however, in contrast to GR, they are not enough in general as they have to be 
complemented with (T2J) and (fit)]) . Hence, the following result has been obtained: 

The proper junction conditions allowing for no shells of matter 
nor branes in F(R) theories with F"(R) ^ are those of GR — the 
agreement of the first and second fundamental forms on both sides of 
the matching hypersurface — together with ([2]) and (flOl) . 

Observe that ([2]) and fflOl) amount to saying that R, as a function, is differentiable 
everywhere (also across E). 

Let me consider the implications of these junction conditions. To start with, 
it should be noted that, in addition to (J32J) one also has ([2]), so that 

n»[R pv } = 0, [R] = (11) 

and R^v has n{n — 1)/2 — 1 allowed independent discontinuities in contrast to the 
n(n — l)/2 of standard GR. One wonders about the allowed discontinuities for 
the energy-momentum tensor. They follow from computing the discontinuity of 
the field equations (IHTj) which, on using (1551) with a = and [K^ v ] = 0, produce 

F'{Rz) [R,A + A h^F'XRv) = k [Tp,] (12) 

where A is a function on E given by (see Appendix) 

A = n^n" [V^V^] 

that is, A represents the discontinuity in the second normal derivative of R across 
E. Taking traces 

(n - l)AF"(ife) = k [T%] . 
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From f llip and ffT2"j) — alternatively from flS]) and 01]) — one derives 

n» [T^\ = (13) 

which happen to be identical with those in the GR case, proving the continuity of 
the n normal components of the energy- momentum tensor across E. Nevertheless, 
in contrast to GR, the n(n — l)/2 allowed independent discontinuities for T^ v are 
not given by those of the Ricci tensor exclusively (as it has one less), but also by 
the new one encoded in the second derivative of R represented by A. In any case, 
f fT3"j) . as well as (J5J) and (TJJ, look like very good properties of F(R) theories. 

2.2.1 Application: perfect fluids 

Consider for example the situation where a perfect fluid interior (V~) is to be 
matched to a vacuum exterior solution {V + ). Then = while 

T^v = QU^Uy + p(g fJiU + u^u u ) (14) 

where g and p are the energy density and pressure of the fluid, while w M is its 
unit velocity vector field (in what follows, I will always use tildes for the energy- 
momentum quantities of the F(R) theories, to distinguish them from the corre- 
sponding GR quantities). Conditions ( 1T3"|) immediately imply then 

pIe = o. ( 15 ) 

£?n^| s = 0. (16) 

The first of this determines the feasible matching hypersurfaces for perfect fluids, 
while the second informs us that the fluid must be either tangent to the matching 
hypersurface or with vanishing energy density there; they are both reminiscent 
of the GR case. 

3 Discussion 

Let me finally discuss the consequences of the above results, in particular in 
relation to whether or not matched solutions in GR are solutions of the extended 
F(R) theories. The answer is generally no, as I am going to argue. I restrict 
myself now to the classical 4-dimensional situation: n = 4. 
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Consider any vacuum solution in GR (with a possible cosmological constant 
A). This is always a vacuum solution of (131]) too provided F(0) = — 2 A. The 
metric is such that 

= A<i R + = AA 

hence 

V M i? + = 0. 

Thus, if it is to be matched to an interior solution, the junction conditions fl2]) 
and (fTUj) imply that, on S, 

iT| E = 4A, V AtJ R"| s = 0. 

This happens to be generically incompatible with an interior metric that matches 
the vacuum solution in GR. The reason is that — letting aside the fact that the 
matching hypersurface might not be the same as in GR — the scalar curvature 
of the interior solution would not be constant on S, even less with vanishing 
derivative, unless in very particular situations. 

Imagine, for instance, that the interior is described by a perfect fluid in GR: 

T ji R = Qu^u y + p(g IM/ + u^Uy) (17) 

where g and p are the energy density and isotropic pressure of the GR perfect 
fluid, respectively. Then, the scalar curvature in the interior region will be given 
by 

R- = -kT gr \ + 4A = k(q - 3p) + 4A 
It follows that only the cases with 

(e-3p)| E = o 

may survive as matched solutions in the extended theories. Notice that then 
£>|s = 0, for the matching hypersurface satisfies p\% = in GR. This already 
discards many global solutions because matched solutions in GR may certainly 
have a jump in the energy density; as a matter of fact, this is what one wishes to 
model in some situations, such as constant density perfect fluids. 

The surviving GR solutions are even more scarce when this is supplemented 
with the last condition 

V>-3p)| E = 

which will certainly rule out many of the few remaining cases. Therefore, most GR 
solutions containing the given vacuum solution as the exterior are not solutions 
of the generalized F(R) theory. 
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4 Explicit examples 



I present two illustrative and important examples. The first one considers the 
general case of static and spherically symmetric stars, the second deals with 
dynamical cases such as collapsing stars, eventually producing black holes. 



4.1 Static spherically symmetric perfect-fluid stars 

Consider any static spherically symmetric line-element 

ds 2 ~ = -A 2 (r)c 2 dt 2 + + rW 

1 — 2m[r)/r 

where r is the area coordinate (round spheres at constant t and r have an area 
of 47rr 2 ) and dVl 2 is the metric on the unit round sphere, m(r) is the so-called 
mass function [24] , and A(r) is a function of r that solves the Einstein's field 
equations for a perfect-fluid energy-momentum tensor (fTTj) with velocity vec- 
tor = —Acd 1 ^. A particular example is given by the interior Schwarzschild 
constant-density solution, see [23] or [22], but there are many others. 

Assume that any such solution has been matched to the exterior Schwarzschild 
solution (I set A = now for simplicity). The matching hypersurface for these 
cases is given by the constant value of the area coordinate r = r such that 

p(r ) = 

and then the mass M of the exterior Schwarzschild part is 

M = — m(r ). 

C Z K 

Eq.([2]) implies then g(r ) = 0, and this by itself forbids many important cases, in- 
cluding the mentioned constant-density solution for the interior region. However, 
even for those cases with a GR equation of state p = p(g) such that g(r ) = 
holds, Eq. ffTOj) still requires that 



; ( r o) = m GR 

dr 

and this of course gets rid of many sensible GR solutions. Actually, apart from 
the case of pure radiation (p = g/3) in GR, only those interiors with 

* W = | W = 0mGR 
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survive. 

Of course, there are many static spherically symmetric interior solutions match- 
ing the Schwarzschild exterior in F(R) theories of gravity. The previous analysis 
just proves that the matched perfect-fluid solutions in GR will not be among 
them generically. 

4.2 Collapsing stars 

As for the second example, I consider dynamical situations. The simplest case is 
given by a Robertson- Walker (RW) interior region 

ds 2 " = -c 2 dt 2 + a 2 (t)dtt 2 k (18) 

where a(t) is the scale factor and dQj, is the complete Riemannian 3-dimensional 
metric of constant curvature k = ±1,0. I denote by g and p the GR energy 
density and pressure of the fluid (|T7|) for the RW geometry, that is to say, [2U [25] 

a 2 + k d 

Kg = 3 ^ — , K(g + 3p) = — 6— 

er a 

where dots denote derivatives with respect to ct. An elementary calculation using 
the field equations ( 131]) proves that the very same metric is also a solution of the 
F(R) field equations for an energy- momentum tensor of a perfect fluid ffl4|) with 
energy density and pressure given by 

Q = + 3p) + ^ + F"(R)3-(g - 3p), (19) 

2 2k a 

p = -P)-^- F"(R){g -3p+-(g-3p)}- F"'(R)(g - 3p) 2 (20) 

2 2k a 

where in these formulae 

/ d 2 + k d\ 

R = K{g-3p) = 6[ 5— + - . 

\ a z a J 

As is well-known, the RW metric (fl8l) matches the Schwarzschild vacuum 
solution in GR whenever p — 0, see e.g. [12], so the fluid becomes dust and 
g(t)a 3 (t) = C =const. The matching hypersurface is co-moving with the RW 
fluid flow. Actually, any such co-moving hypersurface works, and then the exterior 
mass M gets determined accordingly. The Oppenheimer-Snyder collapse (9] to 
form a black hole, see e.g. [251 121], is included here for the closed case k — 1. It is 
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important to remark that the very same matching describes the Einstein-Straus 
model [TT] of a vacuole in an expanding universe, as they are complementary 
matchings in the sense described in [10], see [12J . 

In contradistinction, it is impossible that such RW dust solution matches the 
Schwarzschild vacuum for general F(R). To prove it, observe that the condition 
(F2]) would require 

£?|e = 

but, given that g(t) depends only on t and that £ is timelike, this would lead 
inevitably to 

g(t) = 0. 

Hence, the Oppenheimer-Snyder collapse or the Einstein- Straus vacuole in GR 
are no longer solutions in F(R) gravity. The impossibility of the latter (and other 
RW cavities) in F(R) gravity has been recently obtained in jl]. 

One can further prove that actually no RW space-time can be matched to the 
Schwarzschild solution in F(R) theories with F"(R) ^ (again, I am assuming 
-F(O) = = A for simplicity). To that end, observe that Eq.([2]) implies that 
-R Is = and therefore 

te-3p)| s = o. 

Given that the matching hypersurface is timelike, this implies 

e - 3p = = R~ (21) 

everywhere. Expressions (I19|20p show then that the energy-momentum tensor of 
the RW geometry in the F(R) theory is described by a comoving perfect fluid 
with energy density g and pressure p given by 

g = F'(0)g, p = F'(0)p 

so that the fluid has the same radiation equation of state as in GR: g — 3p = 0. 
But then the matching condition ATS]) implies p(t) = p(t) = 0, and a fortiori 
g(t) = = g(t) too because of (|2~T]) . 

One may wonder if these results depend crucially on the assumption of vacuum 
on the exterior, or on the particular nature of the RW metric. The answer is no 
once again, as similar conclusions follow for general radiating stars with spherical 
symmetry. To describe the exterior of such a star one can use the radiating 
Vaidya metric |T5l [23] 

ds +2 = -(l- ^M] du > _ 2dudr + r 2 dn 2 
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where m(u) is the mass function and u is null retarded time. This is a solution of 
Einstein's field equations for null incoherent radiation (it reduces to Schwarzschild 
for m(u) = M =const.) 





r 



2 dm 
r 2 du 




and it is easy to check that this is also a solution of (13T1) with F(0) = and 



It is known that a very large class of spherically symmetric metrics match the 
Vaidya solution (26] . The only requirement is the existence of a timelike hyper- 
surface £ complying with conditions (|42p . which essentially amounts to finding a 
hypersurface such that n^T~ v is null. The majority of these spherically symmetric 
interiors will have R~ ^ and/or Vi? _ ^ so that they will no longer match 
the Vaidya solutions in F(R) theories. 

Actually, in GR the RW metrics ( 1T81) always match the Vaidya solution and 
the matching hypersurface is timelike (in general not co-moving) whenever the 
dominant energy condition holds: {p/g) 2 < 1 [10J. The mass function gets then 
determined accordingly. The majority of these RW- Vaidya matched models are 
not solutions in F(R) theories, because given that R + = 0, (J2]) implies R~~ = 
and then the argument proceeds as above proving that g = 3p, in which case 



Thus, only the pure radiation RW metric matched to Vaidya is also a global 
solution for arbitrary F(R), keeping the matching hypersurface and the same 
equation of state in the F(R) interior: g = 3p. This is certainly a very meager 
surviving set. Still, it is an explicit example of a matched solution which satisfies 
the GR field equations as well as the field equations of F(R) theories. 
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2 dm 
r 2 du 




g(t)a 4 (t) = const. 
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Appendix 



Let (M ± ,g ± ) be two smooth n- dimensional spacetimes whose respective metrics 
are g . Assume that there are corresponding timelike hypersurfaces E 1 * 1 C 
which bound the regions C on each ±-side to be matched. These two 
hypersurfaces are to be identified in the final glued spacetime, so that they must 
be diffeo morphia The glued manifold is defined as the disjoint union of V + and 
V" with diffeomorphically related points of E + and E~ identified. Henceforth, 
this identified hypersurfce will be denoted simply by E. An indispensable require- 
ment to build a well-defined space-time — with (at least) continuous metric — is 
that the first fundamental forms of E calculated on both sides agree because 
then there is a metric extension g defined on the entire manifold that coincides 
with g ± in the respective V ± and is continuous [27J f!9] . 

In practice, one is given two spacetimes and thus two sets of local coordinates 
{x±} with no relation whatsoever [21]. Hence, one has two parametric expressions 
x± = x± (£ a ) of E, one for each imbedding into each of M ± 1 where {£ a } are 
intrinsic local coordinates for E (/i, u, ■ ■ ■ = 0, 1, . . . , n — 1; a, b, ■ ■ ■ = 1, . . . , n— 1). 
The agreement of the two (±)-first fundamental forms amounts to the equalities 
on E 

Q x ± dx± 

h a b — h ah , h ab = g^ix^)) b . 

Denote by two unit normals to S (one for each side). They are fixed up to a 
sign by the conditions 

±dx± n ± n ±M = i 

and one must choose one of them (say n~) pointing outwards from V~ and the 
other {n't) pointing towards V + . The two bases on the tangent spaces 

are then identified, so that one can drop the ±. The space-time version of the 
now unique first fundamental form is described by the projector to E 

hfii/ Qiii> n^Tii/. 

Notice that 

dx^ dx u 
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At this stage, the Einstein field equations of GR are well-defined in the dis- 
tributional sense, because one can easily prove [HJJ [27J that the Riemann tensor 
distribution (distributions are distinguished by an underline) takes the explicit 
expression 

R a /3iML> = (1 — 9)R a Pnv + 9_R +a p^ v + 5^H a (22) 

where R ±a p^ v are the Riemann tensors of V ± respectively, 9 is the distribution 
associated to the function that equals 1 on V + and vanishes on V~, and 5 s is a 
scalar distribution (a Dirac delta) with support on X acting on any test function 
Y by returning the value of the integral of this function on E: 



It should be observed that 



(f,Y) = [ Y 



V„ 9 = n„, 5 s . 



H a p fll/ in ( l2"2"j) is called the singular part of the Riemann tensor distribution for ob- 
vious reasons, and should only be retained in idealized cases such as braneworlds 
or thin shells of matter. It has an explicit expression 

H a p\^ = -n a [Kfo] n x + n a [Kp X ] - rip [K aX ] % + rip [K afl ] n x (23) 

in terms of the jump of the second fundamental form across sH 

Observe that each second fundamental form is symmetric and orthogonal to n fl , 
thus only the n(n — l)/2 components tangent to S are non-zero. A convenient 
formula for these components is 



From ( |23|) the singular parts of the Ricci tensor Rp and scalar curvature R 
distributions are easily computed to be, respectively 

H p p pll = Hp^ = - [K Plt \ - [K%\ npn^ H% = H = -2 [K%] (24) 



2 Here the standard notation for discontinuities is used, so that for any function / with 
definite limits on both sides of S, one sets for all p € E: [/] (p) = lim x _>j, f + (x) — lim x _>p /~ (x), 
being the restrictions of / to respectively. 
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from where the singular part Qp^ of the Einstein tensor distribution 
G-Pn = Rp P - \g^R = GU + G^(l -6) + 

follows [21] 

g Pli = -[K Pft ] + h Pft [K%], n»g^ = 0. (25) 

In GR, via the Einstein field equation, this provides the singular part r^ iu of the 
energy-momentum tensor distribution 

KT Plt = -[Kp li ]+hp li [K" p ], n^ = (only in GR) (26) 

where k is the gravitational coupling constant: this is known as the Israel formula. 
Observe that only the tangent components Q a b and r a b are non-identically zero. 

As a general result, the Bianchi identity V p R a p pv + V p R a p up + V u R a p pp = 
holds in the distributional sense [19], from where one deduces V^G^ = for the 
Einstein tensor distribution. A standard calculation with distributions then leads 
to 

= = n? [Gfr\ 5 s + (Gp^) 

= <f (nP [Gp p ] + v'Gp, - \n»go°{K% + K; a )^j 

where V denotes the intrinsic covariant derivative within E associated to the first 
fundamental form. This implies the following relations (valid in general and well 
known in GR) 

(*J + K ^ pa = 2nV [Gp p ] = 2nV [Rp p ] - [R], (27) 
V^Gp P = -n p h\ [G pa ] = -n p h% [R pa ] . (28) 

In GR, but not in general F(R) theories, these can be trivially rewritten in terms 
of the energy-momentum tensor T pv and its singular distributional part r pv via 
the Einstein field equations. 

An important remark is that these equations can also be obtained by using 
part of the Gauss and Codazzi equations for £ on both sides, specifically [21] 

R+ - 2R%Tfn v = U~ (K ±p p f + K%K ±P \ (29) 
n»R%h p v = W"K% - V u K ±p p (30) 
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where 1Z is the scalar curvature of the first fundamental form of S. 

Consider now the general case of F(i?)-theories. The field equations read [2] 

F\R)R, V - l -F{R) 9llu - V,V U F'(R) + g pu V p V p F'(R) = kT^ (31) 

where primes denote derivatives with respect to R. Alternatively, these equations 
can be written as follows 

F'{R)R llv - l -F{R)g, u - F"(R) (V P V U R - g pu V p V p R) 
-F"'{R) (V P RV U R - g pu V p RV p R) = kT pv 

where we can see that covariant derivatives of R (up to the second order) are 
needed to compute (13T|) . and the only way to do that for a distribution R with a 
singular part or for a possibly discontinuous function R is by the use of distribu- 
tion theory. One checks that, unless F"(R) = — which corresponds to standard 
GR — , there are terms of type V M V U R and V^RV U R which involve singular terms 
of type 

F(R)g pu , F'(R)H, U , F"(R)V,V U (H5 E ), F" \R)V .([R]^) 
and also 

F'"{R) (v M (#<f ) + [i?K<f ) (v,(#,f ) + [i?K<f ) 

containing products that are not allowed distributional and that can never cancel 
with each other in (|3T|) . In order for the equations (|3T|) to make sense even in 
the distributional sense, a quick analysis implies that the singular part of the 
scalar curvature distribution must vanish H = 0, and from the second in (J24l) 
this entails inevitably 

[K%] = . 

Similarly, the requirement 

[R] = 

is unavoidable unless F"'{R) = 0. There arises an exceptional case for quadratic 
theories where the discontinuity of the scalar curvature can be non-zero and still 
the field equations make sense in the distributional sense. I analyze this case 
separately. 

3 There are some recent results to treat the problem of multiplying distributions, see e.g. 
[28j . but still there is no unequivocal answer to this problem. 
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The generic case F"'{R) ^ 0. 

Taking ([1]) and ([2]) into account, the only remaining singular parts in the field 
equations (l3Tj) are those coming from the Ricci tensor distribution — the first in 
(j24p — and the singular part of 

V U V P R = (1 - 9)V U V^R- + 0V u V p R + + [V^R] n v 8* . 

Taking ([2]) into account, the discontinuity [V^-R] is easily computed [19J to give 

[V P R] = an p , a = n" [V^R] (32) 

where a is a function defined on E that represents the jump in the (normal) 
derivative V-R. Thus, the singular part of the field equations becomes (compare 
with and see 0) 

kt^ = -F'(Rx) [KfJ[ + F"(Rv)n» [V P R] h^, /r p = 

where i?s denotes the value of R at E. 

Using (pQ) and (j2lh the relations ( l2"TH2"Sj) become now 

(K+ + K-JQ" = (K+ + K- a W° = 2nV [R^] , 

= V P H^ = -n p h\ [R pa ] . 

Even though these are different from the GR case, one can still wonder about 
the versions of (12711281) that involve the energy-momentum quantities. To derive 
them, the divergence of (jHJ) is 

-F'(Rx)n p h% [R pa }-F"(Rx) [K pll }V P R^+F"(R^V p a+F"'(R^aV fl Rj: = kTt, 
while from (I3T]) the following discontinuity can be computed 
Kn p h\ [T pa ] = F'(Rx)n p h% [R p<7 ] - F"(R s )n p h% [V P V CTJ R] - F"'(R s )aV p Rx • 
A straightforward standard calculation (19] leads to 

+n v (V M a - [K pp ] V P i? E ) + ^(K+ + K~ u ), (33) 

where A is a function on E defined by 

A = n p n v [V P V U R] . 
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Combining then the previous three expressions one arrives at 

V% M = -n"h\ [T pa \ . 

Computing the total normal discontinuity of ( )3T|) and using ( 1331) and ([1]) one also 

gets 

{K% + K; a )r»° = 2nV [!>„] . 

The exceptional case F"'(R) = 0. 

Assume that F"'(R) = 0, or equivalently that 

F(R) =R-2A + aR 2 (34) 

for some constants A and a. Then a discontinuous .R, with [i2] ^ is allowed in 
principle, from a mathematical point of view. Now, instead of fl32|) one has 

[V^] = an p + V^R) , a = n" [ V M i?] (35) 

and instead of ( 1331 

[V^VuR] = An p n v + n p (V„a - [^] (Vi*) s - A^,V P [i?]) 

+n„ (V M a - (Vi2) E - V*[i2]) + aK% + V„V„[i2] (36) 

where now 

Jfe = \{R + + iT), (Vi2) E = ±V(f2 + + iT)| s , tfJL = + K~ u ). 

The second derivative of the Ricci tensor distribution has then a singular part 
of the form 

v„ ([i?] ) + [v^K-f . 

which, after a calculation with distributions using ( I3"5"j) can be shown to adopt 
the explicitly symmetric form 

5 s | (a - n M n, + + ) + n^ v [R) + n^fl] } + (37) 

where A is a 2-covariant symmetric tensor distribution whose components are 
defined, acting on any test function Y, by 

(A^Y) = - JjR]n, ni/ n^ p Y. 
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It should be observed that this distribution has support on £ but it is of '5" type, 
and thus its product with objects defined exclusively within £ is not defined unless 
extensions of those objects off £ are considered. 

The singular part of the left-hand side of the field equation fl3Tj) reads then 



- {1 + a{R + + RT)}[K^ + 2a% v 
-2aS E | -a V + \ [R] {K+ + K~ v - 2K%n ll n v ) + n^V u [R] + n v V M [R] j (38) 

where Q is a 2-covariant symmetric tensor distribution given by 

so that, acting on any test function Y, one has 

(Q iU ,,Y) = - J[R\h llv n"V p Y. (39) 

From (|38p and the field equations fl3Tl) one observes that the energy-momentum 
tensor distribution is now allowed to have a singular part with several terms 



= T+9 + T~(l -9) + V 6 s + {r^n v + r u n^ + rn^ + 



where 



«r A , = -2aV M [i2], n% = 0, 
«r = 2a [R]K p p . 



Absence of thins shells or branes 

In proper cases where the matching hypersurface intends to describe the boundary 
between a matter interior and a vacuum exterior, or other similar separations, 
only jumps in the matter content (in the energy density, for instance) are to 
be allowed. Thus, one has to require that the singular terms in the energy- 
momentum tensor distributions vanish. From ([3]) or from (JMHJ) this leads to 

[Kfr] = 0- (40) 
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and to 



F"(R s )[R] = 0, 



F"(Rx)nP[V p R} = 0. 



It is easily proven that (|40|) is equivalent to the vanishing of the whole Riemann 
singular part H a ^ lJjU [TP]. If this holds, one can further prove that 

• a local coordinate system around E can be constructed such the metric is C 1 
[T9"] . These are called admissible coordinates [22], but they are not usually 
computed nor, indeed, explicitly constructible in analytical manner. 

• the Riemann tensor has a discontinuity described, in general, by the formula 



where Bpx is a symmetric tensor defined only on E which can be chosen to 
be tangent to E, that is to say, n^Bpx = 0. Thus, Bpx contains n(n — l)/2 
independent components, which are the allowed independent discontinuities 
of the Riemann tensor across E. 

One can easily derive from (14 ip (or alternatively from ( I271|2"B"]) ) that the disconti- 
nuity of the Einstein tensor always satisfies 



so that its n normal components must be continuous across E. Actually, the 
n(n — l)/2 independent discontinuities of the Riemann tensor can always be 
chosen in GR to be those given by the discontinuities of the tangent components 
of the Einstein tensor (and thus, in GR, of the energy- momentum tensor). This 
is different in F(R) of theories with F"(Rj^) ^ 0, as then [R] = necessarily. 
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